Scattering by linear defects in graphene: a continuum approach 
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We study the low-energy electronic transport across periodic extended defects in graphene. In the 
continuum low-energy limit, such defects act as infinitesimally thin stripes separating two regions 
where Dirac Hamiltonian governs the low-energy phenomena. The behavior of these systems is 
defined by the boundary condition imposed by the defect on the massless Dirac fermions. We 
demonstrate how this low-energy boundary condition can be computed from the tight-binding model 
of the defect line. For simplicity we consider defect Unes oriented along the zigzag direction, which 
requires the consideration of only one copy of Dirac equation. Three defect lines of this kind are 
studied and shown to be mappable between them: the pentagon- only, the zz{558) and the 2z(5757) 
defect lines. In addition, in this same limit, we calculate the conductance across such defect lines 
with size L, and find it to be proportional to kpL at low temperatures. 

PACS numbers: 81.05.ue, 72.80.Vp, 78.67.Wj 



I. INTRODUCTION 

Graphene growth by chemical vapor deposition (CVD) 
on metal surface^^M^ is a very promising scalable method 
for producing graphene sheets. However, the present sta- 
tus of the method, typically results in the synthesis of 
polycrystalline graphene abundant in topological defects, 
grain boundaries (GBs) being, by far, the most common 

Due to graphene's hexagonal structure, pairs of pen- 
tagons and heptagons, named Stone- Wales (SW) defects, 
as well as octagons, are expected to form at graph ene 
GBs .El Recent atomic resolution TEM studiePSMSI al- 
lowed the observation of GBs in CVD-grown graphene. 
These experimental studies have shown that the GBs 
are generally not perfect straight lines, and that the 5-7 
defects along the boundaries are not periodic. Further- 
more, as shown by recent TEM studies ,13^ these extended 
pentagon-heptagon defect-lines intercept each other at 
random angles, forming irregular polygons with edges 
showing a stochastic distribution of lengths. This renders 
theoretical studies of such defects difficult, in particular 
when using microscopic tight-binding models. 

Theoretical studies have argued that GBs strongly in- 
ffuence the prop erties of graphene, namely its chemical, 
mechanicaPSEll and electronic ones. Electronic mobili- 
ties of films produced through CVD are l ower t han those 
reported on ex foliated graphene, becausd^^^^ ^ elec tronic 
transporifi^'^ is hindered by grains and GBs.l^^^^ 

Recently, the observation of a linear extended defect 
acting as a one-dimensional conducting charged wird^ 
stimulated some theoretical studies concentrated on the 
scattering and transport properties of such wire.'22H22l 
Most of these studies have so far been focused on tight- 
binding models. 

The use of the continuum approximation on the scope 
of graphene have led to a better understanding of many 
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Figure 1: (Color online) Scheme of a graphene sheet with 
a pentagon-only defect line along the zigzag direction. The 
primitive vectors are ui = a(l, 0) and U2 = a(— 1/2, \/3/2). 



important phenomena occurring in graphene. Moreover, 
we believe that the use of this approach in the study 
of the electronic scattering across extended defects in 
graphene, may further extend our insight onto the physics 
underlying these events in graphene. 

As is widely known, a continuum approximation of 
graphene's first neighbor TB Hamiltonian for states in 
the vicinity of the Dirac points, describes graphene's low 
energy charge carriers as massless Dirac fermions. These 
are governed by two copies of Dirac Hamiltonian, each 
one of them valid around each of the Dirac points. ^'^ In 
this continuum limit, the finite width defect line turns 
out to essentially act as a one-dimensional (infinitesi- 
mally thin) line, separating two distinct regions governed 
by Dirac Hamiltonian. The defect line is modeled by 
a boundary condition on the Dirac spinors, imposing a 
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discontinuity across the defect. This boundary condition 
determines the scattering properties of the defect. 

For simphcity, in this text we only consider extended 
hne defects oriented parallel to the zigzag direction. In 
these cases, we can ignore intervalley scattering, and thus 
consider only one copy of the Dirac equation. While some 
general properties of the boundary condition and trans- 
mittance can be obtained exclusively from the continuum 
description, the specific boundary condition must be de- 
rived from the TB model of the defect. Nevertheless, 
we feel that this approach adds considerably to the un- 
derstandin g of the low energy limit obtained from a TB 
description |22l in particular, it explains, as we will show 
later, why different defects can show exactly the same 
low energy transmittance. 

To illustrate the main physical issues and the method 
of approach, we start with a simplified version of a de- 
fect line, composed of a double line of pentagons oriented 
along the zigzag direction of the graphene lattice (see Fig. 
[ij, which we dub as pentagon-only defect line. Compared 
to the more realistic linear defects that we treat later in 
the paper, the zz(588jP^and the zz(5757) defects, it has 
the added simplicity of full translation symmetry along 
the defect direction, whereas the latter display a dou- 
bling of the unit cell along that same direction. The low- 
energy boundary conditions associated with these defect 
lines are also computed and compared with that of the 
pentagon-only defect line. Suitable choices of the micro- 
scopic parameters lead exactly to the same transmittance 
as a function of angle of incidence in all three cases. Fi- 
nally, we also compute the conductance across a defect 
of length L and find it to be proportional to kpL at low 
temperatures, for all three defects considered. 



II. ELECTRON TRANSPORT ACROSS A 
PENTAGON-ONLY GRAIN BOUNDARY 

A. The continuum description 

A graphene plane with an extended line defect can be 
viewed in the low energy limit as two half-planes of mass- 
less Dirac Fermions, which cannot be joined smoothly, 
because of the defect, a line of discontinuity. To ap- 
proach this problem, consider a finite strip of width W 
in the y direction, where there is a general local poten- 
tial, Viy) = Vs + V^a, + VyCTy + V,a„ for \y\ < 
such that W X (Vs,V) — > (ws,v) as — ^ 0. Integrat- 
ing the Dirac equation in the y coordinate, the resultant 
general boundary condition for the Dirac spinor has the 
form (see Appendix [A| 

*(a;,0+) = A^*(a;,0-), (1) 

where the 2x2 matrix reads 

J^^^-^ay(v,+^.■^)/vp^ (2) 

and <T = {(Tx,c7y,az)- This boundary condition has to 
satisfy the conservation of current in the y direction. 



i.e.,<If'<{x,0+)cry<I>{x,0+) = ^'<{x,0-)ay^{x,0-) for any 
spinor, which implies that A4^ayA4 = ay] the form given 
in Eq. |2] satisfies this condition. An important feature, 
borne out by the derivation of Appendix [X] is energy 
independence of the boundary condition. When we inte- 
grate the Dirac equation across the strip, and take the 
limit W ^ , the term containing the energy e of the 
state, which, unlike the potential, is fixed, drops out. 

An incoming wave from y = —oo, will be partly re- 
flected and partly transmitted at the defect. As a conse- 
quence, the real-space wave- function on each side of the 
defect line is given by 
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where v — ±1 specifies the Dirac cone, 6q is the complex 

phase of lyq^ -\- iqy, and 9^ — — ^q, the complex phase 
of lyq^ — iqy (see Fig. l3|. The sign of the energy is 
noted by s. Imposing the general boundary condition 
gives immediately the following general expression for the 
transmission probability 



T^iE.e) = 



4sin"^ 



(4) 

where we used the property |dct7M| = 1, which follows 
from the condition of fiux conservation, J\4^ayAi — ay. 
A noteworthy feature, that follows naturally from this 
formulation, is the energy independence of the transmis- 
sion across the defect. 

To determine the actual values do the matrix elements 
of M for a specific defect in a graphene lattice we must 
consider its microscopic description. 



B. The low energy limit of tight binding 

The first-neighbor TB Hamiltonian of graphene with a 
pentagon-only defect line (see Fig. [iJ, can be written as 
the sum of three terms, H — + -\- , where 
{H'") stands for the Hamiltonian above (below) the de- 
fect line, while the remaining term, , describes the de- 
fect line itself. In second quantization the explicit forms 
oi and read 



m n ^ 

+ 6^(m— l,n— 1) a{m^n) ^ h.c.y ^ (5) 
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where t is the usual hopping amplitude of pristine 
graphene and is the hopping amplitude between the 
Djn atoms of the defect line, as represented in Fig. [T] 

If we Fourier transform the Hamiltonian along the 
zigzag direction (a;-direction) , we reduce it to an effec- 
tive one-dimensional chain with two atoms per unit cell 
and a localized defect at its center (see Fig. [2]). 




n = —2 n = — 1 n — n = 1 n = 2 
► 

U2 



Figure 2: (Color online) Scheme of the one-dimensional chain obtained by Fourier transformation on the x-direction of the TB 
Hamiltonian of a graphene layer with a pentagon-only defect line along the zigzag direction. The complex hopping amplitude 
t' has the value t' = t{l + e**^""). 



where for {H ) n > 1 {n < —\) . The term describing 
the defect, _ff^, is 

rn ^ 

+ tbHm,0)^d{m) + h.c.\, (6) 



The Hamiltonian of the effective chain is defined as 

H{k^)^H''{k^)+H''{k.,.)+H'^{k^), (7) 

where the three terms on the right hand side of Eq. ([t]) 
read 



^^(fc.) 



E 



[t'bHk,,, 



h.c. 



a{k,j., n 

= —2S,tcos{kxa)d) {kx)d{kx) 
+ th\kx,Q)d{kx) + h.c 



(8a) 

ta){kx,Q)d{kx) 
(8b) 



The one-dimensional chain has alternating hopping 



amplitudes between the atoms, t and t' = til 



Moreover, the electron at a Dm atom acquires an on-site 
energy term, 'e(kx) = —2S^tcos{kxa), which depends on 
the value of the longitudinal momentum k^ ■ 

At the bulk of the one-dimensional chain (n < — 1 and 
n > 0), the TB equations for unit cell n involve ampli- 
tudes at three different positions, rt — 1, n and n + 1, 

eA{kx,n) = -tB{kx,n)-{t'yB(k,,n-l), (9a) 
eB{k^,n) = ^tA{k^,n) -t'A{k^,n+l). (9b) 



Nevertheless, replacing n — ?> n -f 1 in Eq. ( |9a[ ), we can 
solve these equations for A{kx,n+1) and B{kx,n+1) and 
recast them as a recurrence relation relating amplitudes 
at unit cell n + 1 with those at unit cell n, 



L(n + 1) = Tie,k^a).L{n), 



(10) 



where L(n) = [A{kx,n), B{kx,n)]'^ . The passage matrix, 
T(e, 0), is given by 



T(e,0) 



2 cos I 



4 cos 



(11) 



The eigenvectors of matrix T(e, 0) with eigenvalues 
with |Ap = 1, correspond to Bloch solutions propa- 
gating along the one-dimensional chain (band states). 
The eigenvectors with eigenvalues |Ap 7^ 1 correspond 
to evanescent states which decrease when n +00 
(71 ^ -00) when |A|2 < 1 (|A|2 > 1). 

Note that the previous formulation of the TB problem, 
is entirely equivalent to the usual one, where translational 
symmetry along the lattice vectors directions, allows the 
use of Bloch theorem to compute the eigenvectors and 
eigenvalues of the TB Hamiltonian for pristine graphene. 



A similar construction to that of Eq. ( 10 1 can be 



carried out in the rows containing the defect. The TB 
equations for the defect and its neighbors in the one- 
dimensional chain are easily read from Fig. |2] 

e+A{k,,l) = -{t')*B{k,,0)-tB{k,,l), (12a) 

eB{k,,0) = -Wik^)-t'A{k^,l), (12b) 
eD{k,) = -t{A{k,,0) + B{k,,0)) 

- 2(tcos{k^a)D{k^), (12c) 

eA{k^,0) = -it')*B{k,,-l)~tD{k,), (12d) 

e_B(fc„-l) = -t'A{k,,0) -tAik,,-!), (12e) 

where e± = e ± eAV/2, to account for a possible poten- 
tial difference between the two grains separated by the 
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Figure 3: (Color online) (a) Graphene FBZ with the incident 
vectors used in the TB and CA formalism: q = k — K^. (b) 
Scheme of the electron scattering through the barrier (in the 
low-energy limit). 



pentagon- only defect line. This set of TB equations can 
be used to construct a matrix equation relating the TB 
amplitudes in opposite sides of the defect. The technique 
is to solve each equation for the amplitude of the right- 
most site in Fig. [2] and then cast them as 2 x 2 matrix 
equations. For instance, Eq. (12a I is equivalent to 



(*')* 



t t 

1 



Aik^A) 
B{k,,0) 



With this procedure, one can derive 
L(l) = ML(-l), 

where 



(13) 



(14) 



M = i?Afi(e+,0)M2(e,0)M3(e»Mi(e»M2(e_,0)i?^ 

(15) 

is a 2 X 2 matrix; R is the Pauli matrix, used to switch 
rows, (f> = k^a, and A/i, M2 and M3 are 



Mi(e,0) 

M2 (€,</)) 

M3(e,0) 



(1 + e-"t') 




-(l-t-e'"^) 



e+2t£, cos(0) -, 
t 

-1 



(16a) 
(16b) 
(16c) 



Note that the 2x2 boundary condition matrix, M [see 
Eq. (14l], depends on the energy, e, on the longitudi- 
nal momentum, k^, and on the potential difference AV, 
through e_|_ and e_. 

We now have all the ingredients needed to compute 
the scattering coefficients of an electron wave by the 
pentagon-only defect line. Given an incoming wave from 
n = —00, the presence of the defect (line) at rt = 0, pro- 
duces a reffected and a transmitted component. In such a 
case, the wave-functions on each side of the defect (line) 
are given by 



L(n < 0) 
L(n > 0) 



(17a) 
(17b) 



where p and r are, respectively, the reflection and trans- 
mission scattering amplitudes, and 'S'> and '!'< stand for 
the right and left moving eigenstates of matrix T(e, fcj.a), 
the passage matrix for pristine graphene, with corre- 
sponding eigenvalues noted by A> and A<. Imposing 
the boundary condition, Eq. (14), it is straightforward 



to obtain the coefficients p and r for a given energy and 
a given longitudinal momentum. In particular, r reads 



detl 



(18) 



where M = U~^M.U is the boundary condition matrix 
[see Eq. ( 14 1[ in the eigenbasis of the passage matrix of 
pristine graphene T{e,kxa). The transmission probabil- 

2 ~ ^ . 

ity is given by T = |r| — 1/ M22 , since flux conserva- 
tion again requires that |detM| = 1. 

But our main concern is the low energy limit. In the 
following we assume AT^ = 0. Let us consider in parallel 
the equations that propagate the state in the bulk and 
in the defect: 



A(fc^,n-f 1) 
B{k^,n + 1) 

' A{k,,l) 
B{k,,l) 



= T(e,0) 
= M(e,</.) 



A{ka;,n) 
B{kx,n) 

A{k,,-l) 
B{k^,~l) 



bulk; (19a) 
defect; (19b) 



As is well known, near a Dirac point K^, the slowly vary- 
ing Dirac spinor 5''^(r) is deflned by (ignoring irrelevant 
normalization constants) 



1'''(mui +nU2) = e-'^''-('""i+""^ 
and for a plane wave along Ui 
^""(toui -I- nu2) 



A(m, n) 
B(m, n) 



(20) 



A{kx,n) 
B{kx,n) 



(21) 



where q = k — K^. This allows us to recast Eqs. ( |19a l 
and ( |19b[ ) in terms of the Dirac flelds, 

(g„ (n + 1)U2) = e-*K'' "^T(e, 0)*''(g,, nua), (22a) 
*'^(fe, U2) = e-*K-2u.j^(^^ 0)*''(fe, -U2), 

(22b) 

where • U2 = — i^27r/3. 

If we take the Fourier transform with respect to the 



spatial variable along U2 in Eq. ( 22a I , 



(23) 



In Appendix |B] we show that the matrix multiplying 
on the right hand side tends to the identity matrix when 
q, e — 7- 0; if we expand the right hand side to linear 
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order in e and q, we obtain, as we should, the Dirac- 
Weyl equation (see Appendix [B|. However, at the defect, 
we find 



e (rad) 



e (rad) 



-iu27r/3v 





-1 



when (f), e ^-0, (24) 



which gives rise to the following equation 



*"(9.,U2) 



1 

-1 e 



*"(9.,-U2). 



(25) 



After Fourier transforming the previous equation in 
Qx, and as the continuum approximation yields a — in 
U2, near the Dirac point, we end up concluding that the 
defect introduces a discontinuity in the Dirac fields of the 
form we derived from general considerations, '^'^{x, 0^) = 
7W«'''(a;,0-), with 



M 



1 

-1 e 



(26) 



The transmission probability, given by the general ex- 
pression of Eq. (|4| , becomes here 



(27) 



In Fig. [4] we plot the transmission probability T, in 
terms of the angle of incidence on the defect line, for both 
the TB and the continuum approximation (CA), with 
^ = 1.2. The various plots refer to different energies, but 
always to the same Dirac point {v — 1). As expected, 
the lower the energy, the better the agreement between 
the TB and the CA results. For the other Dirac point, 
the results are mirror-symmetric relatively to the normal 
incidence angle 6 = tt/2. 

A special case is of some interest, namely, for low en- 
ergies, AV = 0, and ^ = 2, the transmission probabil- 
ity becomes T^{9) = (1 -I- vcos9)/2, in which case the 
pentagon- only defect line acts as a valley filter, for angles 
of incidence close to 6' = 0, vr. This same feature has been 
found in another type of defect, the zz(558), which we 
consider in the next section, by Gunlycke and WhitCj'^il 
this is no accident; we will show that these two defects 
share the same low energy limit. 

It is worth noting, that since the passage matrix in 
the continuum limit is obtained with (j) = Vi^ ■ Ui and 
e = 0, it is easily got in a hack of envelope calculation, 
by writing and solving the TB equations at zero energy. 
This procedure is carried out in Appendix [Cj 

It is expected that defect lines and grain boundaries 
in graphene are reactive,!^ being a likely location for ad- 
sorption of atoms or molecules. Such adsorbates, are 
expected to locally perturb the properties of the defect 
lines. For simplicity, we may assume that the adsorbate 
only modifies the local energy at the atom it adsorbs to. 
We can account for such a phenomenon in the pentagon- 
only defect line, including in its TB model, an on-site 
energy, eq at the D atoms of the defect line (see Fig. [l] 




I (rad) 



(rad) 



Figure 4: (Color online) Plot of the transmission probabihty, 
T, in terms of the angle of incidence on the pentagon- only 
defect line, 6^ (for the low-energy limit, around K^, with 
V = -1-1). We have used the value ^ = 1.2 to obtain these 
curves. In each of the panels, we compare the TB result for 
e > (full violet curves) and for e < (dashed violet curves) , 
with that obtained from the CA (in green), with AV^ = 0. 
Panel (a), (b), (c), (d), (e) and (f), stand, respectively, for 
energies \e/t\ = 0.01, \e/t\ = 0.04, \e/t\ = 0.08, \e/t\ = 0.16, 
\e/t\ = 0.32 and \e/t\ = 0.64. 



or Fig. |2]). Such a modification of the TB model, will 
necessarily modify the TB boundary condition matrix. 



M [see Eq. (14)], as well as the continuum approxima- 
tion one, [see Eq. (26l]. The TB boundary condition 
matrix, [see Eqs. ([T6|], will have its e + 2£tcos(kxa) term 
modified. This will now include the on-site energy, eq, as 
e' — e-\- 2^t cos{kxa) -\-eo. In the CA limit, the boundary 
condition matrix, Ai, will have {S^t — eo)/t in the A^22 
entry of the matrix instead of ^. Thus, the adsorption of 
molecules at the defect line, in very low energies, will be 
equivalent to rescaling the hopping between the D atoms 
at the defect line. 



III. THE 2z(558) AND THE zz{5757) DEFECT 
LINES 

We now extend this treatment to the case of a zz(558) 
defect linei"2i22 (g^g ^ig. [s]), and of a zz(5757) defect 
line (see Fig. |6|. 

We can proceed in close analogy with the case of a 
pentagon-only defect line treated in the previous section. 
But these more realistic defects exhibit a feature that is 
not present in the previous case, namely, the doubling 
of the unit cell in the direction parallel to the defect. 
The corresponding folded First Brillouin Zone (FBZ) has 
twice as many states at the same Bloch wave vector, as in 
the original FBZ of graphene; the real space unit cell has 
two A {Ai, A2) and two B {Bi, B2) sites. Around the 
new Dirac points, now located at K± — ±7r/3(l, — -\/3), 
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m - 1 rn m + 1 
Figure 5: (Color online) Scheme of a Z2:(558) defect line. 




- 1 m m + 1 

Figure 6: (Color online) Scheme of a zz(5757) defect line. 



there will be, in addition to two low-energy Dirac cones, 
two high energy bands.'^^l At low energies, e ~ 0, the 
extra states show up as evanescent solutions .1^^^ 

In pristine graphene we know the form of the high and 
low energy modes since they are Bloch states of different 
wave vectors in the unfolded Brillouin zone. We can use 
this to define a change of basis that decouples, in the 
bulk, these two energy sectors (0 = kxo)'- 







[All 
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= m 




(28) 
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(29) 



Defining 

L(n) = [A+{k.^,n),B+{kx,n),A^{kx,n),B^{kx,n)]'^, 



(30) 



we have, L(?i + 1) = Tc;L(?i), where the matrix Td, writ- 
ten in blocks of 2 x 2 matrices, is 



Td(e,0) 



T+(e, 






T_(e,, 



(31) 



The -|- and — amplitudes propagate independently; 
and T_, the passage matrices associated with the high 
and the low-energy TB modes, are 



T+(e,(/.) 



2cos(f) 



2zsin(|) 



1 



4cos2(|) 



4sin^(|) 



(2 



(32a) 
(32b) 



The above computations are due to appear in a compan- 
ion papei^^ devoted to the study of these same systems 
under the TB approach. 

In parallel with what we have done for the pentagon- 
only defect line [see Eqs. (12 1- (16l], using the TB equa- 
tions at the zz(558) or at the zz(5757) defect lines, it is 
possible to write an expression relating amplitudes at the 
two sides of these defects, L(l) = ML(— 1). The matrix 
M is now a 4 X 4 matrix relating the four amplitudes 
at each side of the defect line, and admixing, in general, 
high and low-energy modes of different sides of the de- 
fect. The high energy sector passage matrix in the bulk 
near e = and (p = • Ui = vtt/S is 



— _„-i7r/6 







1 

x/3 



(33) 



The corresponding eigenstates are evanescent, one grow- 
ing exponentially as e'-"'°S'^)/^, localized on the B sub- 
lattice, and the other decreasing as e^*^"'°S'^^/^, localized 
in the A sub-lattice. This same result was obtained by 
Ostaay et al. in the scope of the total reconstruction of 
the zigzag edge by Stone- Wales defects .^^Sl 

Given this, we conclude that a low energy state, must 
have the following form in each one of the sides of the 
defect 



^{k^,n) 





B+{k^,n) 
A_{kx,n) 
B^{k^,n) 

A+{kx,n) 


A-{k.^,n) 
B-{k^,n) 



n<Q; (34a) 



n> Q. 



(34b) 



This form fixes the i3+(fca;,— 1) amplitude, in terms of 
the low energy amplitudes ^_(fca:,— 1) and _B_(fca;,— 1), 
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since 



1.0 



ll22B+(fc„-l)+M23A_(A:„-l)+M24B_(A:„-l) = 0, qS 

(35) 



and leads to an effective boundary condition relation for 
the low energy amplitudes only. The latter reads 



(36) 



where the effective boundary condition matrix is ob- 
tained from matrix M 







■ A_(fc„-1)" 


B_(fc„l) 




S_(fc,,-1) 



ircff 



n33-M32M23/M22 M34 
143 -M42M23/M22 M44 



II32M24/M22 
II42M24/M22 



(37) 



The low energy sector, with the matrix T_(e, 0), can be 
analyzed exactly as was done in Appendix [B] for the pen- 
tagon only boundary. We define the Dirac fields as be- 
fore, 



*"(?:,, nU2) 



-iKiy ■nu2 



A_{ka:,n) 



(38) 



so that, in the bulk 

^%q,, {n + 1)U2) - e-*K„.u.Tp_ 0)vl/-(g,, nU2) 



(39) 



With a procedure entirely similar to the one detailed in 
Appendix [B] one finds, after Fourier transforming in n, 
that ^I^q satisfies the Dirac equation. At the defect. 



vl/''(g„ U2) = e-^K-2"^M°« [Q,v-] vl'''(fe, -U2) 



(40) 



The calculation of yields 



1 
-1 2|| 



(41) 



and so the boundary condition for the Dirac fields is 





-1 



9 £,2 

^1? 



(42) 



It is remarkable that this has exactly the same form as 
found in the pentagon-only boundary (c./. Eq. |26| ); the 
low energy transmission probabilities of these two line 
defects are the same provided 2^2/^1 = ^- In Fig- |7]we 
compare the transmission probabilities calculated with a 
full TB calculation for different values of ^1 and ^2 but the 
same value of ^2/^11'^^ ^^nd the corresponding low energy 
approximation. 
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Figure 7: (Color online) The transmission probabilities for a 
2^(558) defect, at e/t — .03, with ^1 = 1, ^2 = 0.5 (dashed 
red), ^1 = 1.5, ^2 = 1.125 (dashed-dot, blue), obtained in a 
full TB calculation, and the corresponding low energy approx- 
imation (continuous black), given by Eq. (27 1 with ^ = 1.0. 



The treatment of the zz{5757) line defect presents no 
further novelty. Using the quick derivation method out- 
lined in Appendix [C) we arrive at the following passage 
matrix for the Dirac fields 



M 



off 
5757 



-1 



26(4' 



a b 

-b c 



(43) 



where a = - H/^), b = -U^b " Ca) and c = 

2(^b + d + iUD/ic (see Fig [g] for the notation of the 
hopping amplitudes). The form of the passage matrix 
(and the transmission probability) is not identical to the 
previous cases, unless — — ^q/2. 

Despite the similarities, for general values of the hop- 
ping parameters, the transmittances originating from 
each of the defect lines can be considerably different. 
Such a case can be seen in Fig. [8] where we compare 
the low-energy transmission probabilities associated with 
each one of the previously discussed defect lines, for a 
special situation with all hopping amplitudes equal to t, 
the bulk nearest neighbor amplitude. 



IV. CONDUCTANCE 

In this section we address the calculation of the linear 
conductance, across a line defect, and show that the en- 
ergy independence of the transmission probability at low 
energies found in the previous three cases gives rise to a 
conductance linear in fc^- 

We make the usual assumption that the electron reser- 
voirs at each side of the defect line are in equilibrium 
and are thus described by the single particle Fermi-Dirac 
distribution. Then the expression for the total net cur- 
rent across the defect line, associated with electrons liv- 
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where C" = Lg^ggC^ / {Att^1t?vf)- The transmission prob- 
ability, T{E,/W = 0) = T(0,0), does not depend on E, 
as was seen above; it does depend on the values of the 
hopping amplitudes in the vicinity of the defect, through 
the passage matrix, [see Fig. [9] , for the case of the 
pentagon only defect]. We obtain. 



T(0,0) 



ksT X h 



(47) 

where the function h{x) is a Fermi integral, h(x) := 
I-^ dy \y + x\ exp(y)/ [exp{y) + if, with limits h{0) = 
21n(2), and h{x) — > |x| for ^ 1. The conductance 
is linear in temperature for <C fc^T; in the opposite 
limit, ^ ksT, it is practically temperature indepen- 
dent, 



Figure 8: (Color online) Comparison between the transmis- 
sion probabilities, T, across the defect for: the pentagon-only 
[or Z2(55)[ defect line (full orange line); the zz{558) defect line 
(dashed red line); the zz{5757) defect line (dot-dashed green 
line); The transmission probabilities were calculated in the 
low energy limit from Eq. Q, with all the hopping parameters 
were chosen equal to 1, ^ = = ^2 = = ^6 = = 1- T is 
plotted in terms of the angle of incidence, 6'q , of q = k — K+ , 
the momentum around the Dirac point K-|-. The transmission 
probabilities valid for the vicinity of the other Dirac point, 
K_, are obtained from the ones plotted above, by a reflection 
along the vertical line 6 = n/2. 



ing around the Dirac point K^, is given by 



r^iAv) = c 



/•OO 




/ dE 




' —00 


2 



eAV 

f{E,fi+—)-fiE,,, 



T''{E,AV) 
eAV^ 



(44) 



where C = gsc/ {'iii'^h^vp), gg stands for the spin degen- 
eracy, AV is the potential difference between the each 
side of the defect line (see Fig. [T|), and f{E,fj,) the 
Fermi-Dirac distribution function for chemical potential 
fi: Afi = eAV is the difference between chemical poten- 
tials at the two grains. In addition, T'^{E,AV) stands 
for an angle-integrated transmission probability. 



r{E,AV) 



T''{E,AV,0)smed9, (45) 



The total current is obtained summing the currents asso- 
ciated with the two Dirac points Jy — + Jy . One can 
verify that T+{E,AV) = T-{E,AV) = T{E,AV), and 
thus, Jy — gvJy ; where g^ — 2 is the valley degeneracy. 

The conductance is then given by G = LJy/ AV , where 
L is the length of the defect line, when the current is in 
the linear regime, 

G(T) = C \E\T{E,Q) (- ^^g^^^ ) dg, (46) 



G(T = 0) = Lg.gs 



-T(0,0)|Ai| 



^T(0,0)fc.L. 



(48) 



0.75 




■ 0.25- 



Figure 9: (Color online) Dependence of the total angular in- 
tegrated transmission probability associated with one Dirac 
point, T(S,Ay) ^T+{E,AV) = T-(S, AV), when AV = 0, 
in terms of the hopping, ^t, between the atoms D at the de- 
fect line (see Fig. [TJ. The vertical dashed line, indicates the 
value of ^ used to obtain the curves of Fig. |4] ^ = 1.2. 



V. CONCLUSION 

In this text we have focused on the study of the low- 
energy continuum limit behavior of the electronic trans- 
port across periodic defect lines oriented along the zigzag 
direction of graphene. We have argued that in this limit, 
such extended defects essentially act as one-dimensional 
infinitesimally thin lines, that separate two regions gov- 
erned by the Dirac Hamiltonian. In the low-energy con- 
tinuum limit the defective line imposes a boundary con- 
dition on the Dirac spinors at each side of the defect. 
It is this boundary condition that defines the low-energy 
behavior of the electronic transport of such systems. We 
have demonstrated how can the boundary condition valid 
in the low-energy limit be computed from the TB descrip- 
tion of the defect line. 



9 



We have presented such a reasoning while working out 
the problem of the electronic transport across a pentagon- 
only defect line, finding its transmittance to be energy in- 
dependent. Furthermore, we have also studied two other 
kinds of more realistic periodic defect lines: the zz(558) 
defect line, recently observed in graphene sheets,^ and 
the zz{5757) defect line. We have briefly examined these 
latter cases, emphasizing the fact that their periodicity 
forces the appearance of high-energy modes at the Dirac 
points in addition to the typical low-energy Dirac modes. 
It has been shown that the influence of the former, can 
be encompassed in an effective boundary condition seen 
by the low-energy massless Dirac fermions. The trans- 
mittance originating from such boundary conditions was 
again found to be energy independent. Furthermore, we 
have pointed out that the effective boundary conditions 
arising from the zz(558) and zz(5757) defect lines, turn 
out to be similar to the one arising from the pentagon- 
only defect line. Moreover, the former can be mapped 
into the latter by an appropriate choice of the hopping 
parameters at the defects. 

It is important to note that by expressing the trans- 
mission probability in terms of the boundary conditions 
satisfled by the Dirac flelds at the defect, these results 
cast some light on the low-energy limit of the full TB 
calculations,'^^ leaving us with a better understanding of 
the physics underlying such systems. 

We have in addition shown how can we compute the 
low-energy limit conductance expression, across this kind 
of defect lines. Interestingly, at low temperatures, the 
conductance across a defect line of size L, turned out 
to be linear in kpL. This feature originates from the 
energy independence of the low-energy transmittance of 
our defect lines. 

Finally, we must mention, that the procedures pre- 
sented in this text, can be used to solve more complex 
scattering problems in graphene. Not only linear and pe- 
riodic defect lines, but also periodic curvilinear extended 
defects oriented along graphene 's zigzag direction, can be 
worked out using the framework presented in this text. 
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Appendix A: The general low-energy boundary 
condition for a defect line oriented along graphene's 
zigzag direction 

In this appendix we derive Eq. ([2]), determining the 
general form of the boundary condition matrix of a zigzag 
oriented defect line in the continuum limit. 

Suppose that our defect line is located in the region 
deflned hy y G [0, W]. Assume then, that in this region, 
we have a constant potential term in the Dirac equation 
of the form 



V 



(Al) 



Our aim is to consider the limit where W{Vs,~V) — ?> 
(mo, vq) when — > so that we can obtain the boundary 
condition of the defect line in the continuum limit. We 
must refer that there are some works published on the 
literature, considering the electronic scattering across re- 
gions with p otentials that are a particular form of that 
given in Eq. |AiP^^' 

The Dirac equation in the region of the potential is 



{va^{-id^) -f (Jy{-idy)) + (F, + V • cr) 



(A2) 



Since the defect line is oriented along the a;-direction, we 
can choose 

<^,{x,y) = (A3) 
which, after substitution into Eq. ( |A2| , results in 

dy^^iv) = tP^.iy), (A4) 
where the operator P reads 

P = ^[e-i^VFqx<Jx-V,-V -a] (A5) 

VF 



This flrst order differential equation (A4| can be straight- 
forwardly integrated. 



(A6) 



Taking now the Hmit W{Vb,Y) — >• (mo,vo) when W — )■ 
0, we obtain the following expression for the continuum 
limit of the boundary condition of a zigzag oriented defect 
line 



$.(0+) 



-io-j, (llo +vo • cr) /t)F (|Q- 



(AT) 



just as in Eqs. ([T| and 

As a final comment, we must stress the fact that the 
remaining terms in P, namely e and vpqxO'x, do not 
contribute to the boundary condition when we take this 
limit; they are fixed in value, unlike the potential terms, 
and cannot give rise to a discontinuity when W ^ 0. 
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Appendix B: Dirac equation from the passage 
matrix 

In this appendix, we show how the Dirac- Weyl equa- 
tion can be obtained from the low-energy passage matrix 
relation in Eq. ( 23 1 



(Bl) 



Since 4> = (K^ -I- q) • Ui = vAir/'i + q • Ui, 



47r 



As we are working out a theory valid around the Dirac 
points, q and e are small, and we can thus expand the 
exponential, keeping solely the first order terms in the 
momentum. 



e-"i"^T £,1/ 



47r 



q Ui 



-iq-(ui/2+U2) 



2 cos {y 



2tt_ 
3 



q ui 

2 , 



" f 




-1 0. 





1 4cos2(.f -f^)- 

-iqy - vq^ 

-iqy + vqx 



(B3) 



where / is the 2x2 identity matrix. This term cancels 
the one in the right hand side of Eq. ( B2 1 and we are left 
with 



where t' = <(l + e «'=^«/2), From now on, we set ourselves 
at kj.a — vAii/'i. The five equations written in Eqs. (CI I, 







t 

f 



q 



x/3 



-iqy - vq^ 

-iqy + vq^ 



*q = 0. 



Upon multiplying by i<Ty, we obtain Dirac's equation, 

e^-^ = Wi.(T, -q*^, (B5) 

where cr^ — {vcTx, (^y), with the usual notation of = ±1 
identifying the Dirac point. 



contain 7 amplitudes, and we can solve them all in terms 
of A{kxT—l) and B{kxT—^)\ we obtain for yl(fca;, 1) and 
S(fc.,l), 



(B4) 



' A{kx,l) ' 




1' 




' A{kx,-l)' 


B{kx,l) 




.-1 ^. 




B{kx,-l) 



,(C2) 



Appendix C: Quick derivation of the continuum 
limit boundary condition for the pentagon-only 
defect 



In this brief appendix, we will present a quick deriva- 
tion of the pentagon-only defect low-energy boundary 
condition (e = and q^ = 0) , Eq. ( 26 1 . 



Let us start from the TB equations at the pentagon- 
only defect, Eqs. (12 1. In these equations, we begin by 



setting e = 0. In this way, the TB equations at the defect 
now read 



= 


-{t')*B{kx,0)-tB{kxA), 


(Cla) 


= 


-tD{kx)-t'A{kx,l), 


(Clb) 


= 


-t{A{kx,Q) + B{kx,Q)) 






~2£t cos{kxa)D{kx), 


(Clc) 


= 


-{tyB{kx,-i)-tD(kx), 


(Cld) 


= 


-t'A{kx,0)-tAikx,-l), 


(Cle) 



which, using Eq. (22b I, immediately identifies the pas- 



sage matrix for the Dirac fields, Eq. (26 1 



This procedure is quite general, and can be applied to 
the the other line defects considered in this paper. In 
that case, however, we must express the TB amplitudes 
in terms of the amplitude of the low and high energy 
modes and set to zero the evanescent amplitudes of the 
states that grow on each side of the defect. We can then 
solve for the low energy amplitudes on one side of the 
defect, and obtain directly the 2x2 passage matrix for 
the propagating modes. 



^ X. Li, W. Cai, J. An, S. Kim, J. Nah, D. Yang, R. Finer, 1312 (2009). 

A. Velamakanni, I. Jung, E. Tutuc, et al.. Science 324, 



11 



^ A. Rcina, X. Jia, J. Ho, D. Nezich, H. Son, V. Bulovic, 

M. S. Drcssclhaus, and J. Kong, Nano Lett. 9, 30 (2009). 
^ K. S. Kim, Y. Zhao, H. Jang, S. Y. Lee, J. M. Kim, K. S. 

Kim, J.-H. Ahn, P. Kim, J.-Y. Choi, and B. H. Hong, 

Nature 457, 706 (2009). 
* S. Bae, H. Kim, Y. Lee, X. X. andJae Sung Park, Y. Zheng, 

J. B. andTian Lei, H. R. Kim, Y. L Song, Y.-J. Kim, K. S. 

Kim, et aL, Nature Nanotechnology 5, 574 (2010). 
^ P. Y. Huang, C. S. Ruiz-Vargas, A. M. van der Zande, 

W. S. Whitney, M. P. Levendorf, J. W. Kevek, S. Garg, 

J. S. Alden, C. J. Hustedt, Y. Zhu, et aL, Nature 469, 389 

(2011). 

® K. Kim, Z. Lee, W. Regan, C. Kisielowski, M. F. Crommie, 

and A. Zettl, ACS Nano 5, 2142 (2011). 
^ P. Nemes-Inczc, K. J. Yoo, L. Tapaszto, G. Dobrik, 
J. Labar, Z. E. Horvath, C. Hwang, and L. P. Biro, Appl. 
Phys. Lett. 99, 023104 (2011). 
® A. Stone and D. Wales, Chem. Phys. Lett. 128, 501 (1986). 
® J. Meyer, C. Kisielowski, R. Emi, M. Rossell, M. Crommie, 

and A. Zettl, Nano. Lett. 8, 3582 (2008). 
^° J. Lahiri, Y. Lin, P. Bozkurt, I. I. Oleynik, and M. Batzill, 

Nature Nanotechnology 5, 326 (2010). 
" A. Ferreira, X. Xu, C. Tan, S. Bae, N. Peres, B. Hong, 
B. Ozyilmaz, and A. Castro Neto, EPL 94, 28003 (2011). 
^■^ S. Malola, H. Hakkinen, and P. Koskinen, Phys. Rev. B 
81, 165447 (2010). 

O. V. Yazyev and S. G. Louie, Phys. Rev. B 81, 195420 
(2010). 

1* R. Grantab, V. B. Shenoy, and R. S. RuofF, Science 330, 
946 (2011). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 



Y. Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. 

Firsov, Science 306, 666 (2004). 
^® K. S. Novoselov, D. Jiang, T. Booth, V. V. Khotkevich, 

S. M. Morozov, and A. K. Geim, Proc. Natl. Acad. Sci. 

102, 10451 (2005). 
^'^ A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 

(2009) . 

N. M. R. Peres, Rev. Mod. Phys. 82, 2673 (2010). 

O. V. Yazyev and S. G. Louie, Nature Materials 9, 806 

(2010) . 

^° D. A. Bahamon, A. L. C. Pereira, and P. A. Schulz, Phys. 
Rev. B 83, 155436 (2011). 

D. Gunlycke and C. T. White, Phys. Rev. Lett. 106, 
136806 (2011). 

^2 L. Jiang, X. Lv, and Y. Zheng, Phys. Lett. A 376, 136 

(2011) . 

G. W. Semenoff, Phys. Rev. Lett. 53, 2449 (1984). 
J. N. B. Rodrigues, P. A. D. Gongalves, N. F. G. Rodrigues, 
R. M. Ribeiro, J. M. B. Lopes dos Santos, and N. M. R. 
Peres, Phys. Rev. B 84, 155435 (2011). 
J. A. M. van Ostaay, A. R. Akhmerov, C. W. J. Beenakker, 
and M. Wimmer, Phys. Rev. B 84, 195434 (2011). 
^® J. N. B. Rodrigues, N. M. R. Peres, and J. M. B. Lopes dos 
Santos, to be submitted (2012). 

M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. 
Phys. 2, 620 (2006). 
^® J. Viana Gomes and N. M. R. Peres, Journal of Physics: 
Condensed Matter 20, 325221 (2008). 



